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Abstract 

We discuss a general strategy to compute the coefficients of the QCD chiral La- 
grangian using lattice QCD with Wilson fermions. This procedure requires the 
introduction of a lattice chiral Lagrangian as an intermediate step in the calcu- 
lation. The QCD chiral Lagrangian is then obtained by expanding the lattice 
effective theory in increasing powers of the lattice spacing and the external mo- 
menta. In order to investigate the consequences of the chiral symmetry breaking 
induced by the Wilson term, we study the lattice chiral Lagrangian at the leading 
order of the strong coupling and large N expansion. We show that the effects of 
the Wilson term can be conveniently taken into account, in the lattice effective 
theory, by a suitable renormalization procedure. In particular, we show that, at 
the leading order of the strong coupling and large expansion, the chiral sym- 
metry is exactly recovered on the lattice provided that the bare quark mass and 
the lattice operators are properly renormalized. 
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1 Introduction 



The QCD Chiral Lagrangian, originally introduced by Weinberg |T[| as a convenient 
way to reproduce the predictions of PCAC and current algebra, has been subsequently 
promoted to the role of a consistently renormalizable effective theory p|, ^, which 
provides a powerful tool for describing the phenomenology of low-energy QCD. A recent 
review of the phenomenological results obtained in this framework can be found in 
refs. i-i. 

The basic observation of the chiral Lagrangian approach is that, at very low energy, 
the light pseudoscalar mesons are the only hadronic particles which can be close to the 
mass shell in the strong interaction processes, and can therefore contribute poles to the 
analytic structure of the amphtudes. Once these contributions are accounted for by 
the introduction of a suitable set of fields, an expansion in powers of momenta becomes 
possible. The QCD chiral Lagrangian is then the most general effective theory, con- 
structed in terms of pseudoscalar fields and a set of external sources, which is consistent 
with QCD chiral symmetry and its assumed spontaneous dynamical breaking. 

Another important feature of low-energy strong interactions is that, in the limit 
of massless light quarks {u, d, s), chiral symmetry prevents the interactions among 
the pseudoscalar Goldstone bosons at zero external momenta. Therefore, the effective 
theory can be systematically organized as an expansion in increasing powers of the 
quark masses and external momenta: 

= £2 + A + • • • (1) 

The most general form of the chiral Lagrangian is constrained by the symmetries of 
the fundamental QCD theory. This fixes completely the form of the interactions among 
the pseudo-Goldstone bosons and the correlation functions of the external currents in 
the low-energy limit. However, a set of numerical coefficients, which define the strength 
of these couplings, cannot be fixed by symmetry requirements alone. 

At the lowest order p^, the effective Lagrangian £2 can be expressed in terms of 
only two arbitrary parameters which represent, in the chiral limit, the pseudoscalar 
decay constant F.,^ and the quark condensate. One finds then that ten additional chiral 
couplings are necessary to describe the low-energy QCD phenomenology at 0{p^) 
They are usually denoted by the symbols Li, . . . , Liq. The main purpose of this paper 
is to discuss a possible first principle lattice calculation of these coefficients. 
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At present, our knowledge of the coefficients of QCD chiral Lagrangian mainly comes 
from the experimental data. The couplings are fixed by comparing the predictions of the 
effective theory with a set of physical amplitudes measured in the experiments. Further 
constraints on the coefficients can be derived by general theoretical considerations, 
relying, for instance, on the large N expansion. Such a determination is typically 
affected by large uncertainties, and some of these coefficients are still known with a 
relative error larger than 100% 0. 

From a theoretical point of view, the coefficients of chiral Lagrangian ought to 
be calculable from the QCD Lagrangian. They are functions of the fundamental QCD 
scale, Aqcd, and the heavy quark masses. A rough estimate of their size can be obtained 
by requiring that the tree-level amplitudes from £4 are of the same order of magnitude of 
the loop corrections from C2 H]- This leads to the prediction Li ~ (F^/4)/A^ ~ 2x 10~^, 
where A^^ ~ AttFj^ is the scale of chiral symmetry breaking. Such an estimate reproduce 
in fact the correct order of magnitude of the chiral coefficients renormalized at the 
rho mass scale, indicating a good convergence of the momentum expansion below the 
resonance region. 

A more quantitative picture can be obtained by exploiting the role of meson reso- 
nances in the low-energy effective theory 0. The basic idea consists in writing down 
the most general low-energy effective Lagrangian for effective vector, axial, scalar, pseu- 
doscalar and tensor fields. The standard chiral Lagrangian is then obtained by inte- 
grating out the mesonic degrees of freedom but for the pseudoscalar ones. In this way, 
one ends up with an explicit expression for the chiral coefficients in terms of low-energy 
mesonic parameters, masses and decay constants. This approach provides values of the 
chiral couplings in good agreement with the experimental determinations, showing that 
these couplings are almost completely saturated by the mesonic resonances contribu- 
tions. However, in order to determine their actual values, one must still rely on the 
experimental measurements. 

In principle, the QCD chiral Lagrangian can be obtained by integrating out the non 
effective quark and gluonic degrees of freedom from the original fundamental theory. 
Indeed, being the two theories mathematically equivalent, the partition function Z can 
be expressed as: 

Z = (^dAdipdipj exp < i / d'^x Cqcd \ — (dU) exp < i / d'^x C^jj \ (2) 
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However, although the integral over the fermionic fields in the above equation can be 
performed explicitly, we do not know, in the general case, how to perform analytically 
the remaining integration over the gluonic fields. 

A more feasible theoretical approach would consist in performing a matching be- 
tween the effective and the fundamental theory. Specifically, one could compute a 
sufficient number of physical amplitudes, both in the effective and the original theory, 
and derive the values of the chiral coefficients from a comparison of the results. The 
calculation in the full theory, being non perturbative, would then require the imple- 
mentation of a numerical lattice simulation. 

However, such an approach would be affected by the following technical difficulty: 
in order to perform the matching between the effective and the full theory, one should 
consider a set physical amplitudes defined in the region of low external momenta, p < 
rrip, where the predictions of the chiral Lagrangian can be reliably obtained. On the 
other hand, in the numerical calculation, an intrinsic infra-red cut-off is introduced by 
the finite size of the lattice. The minimum value of momentum that can be considered 
is Pmin = 2-71 /La, whcrc L is the lattice size and a the lattice spacing. In current lattice 
calculations, Pmin is typically of the order of the rho mass and thus lies at the border of 
the kinematical region accessible to the calculations in the effective theory. In order to 
overcome this difficulty, one should consider either larger values of the lattice spacing, 
thus increasing the finite cut-off effects on the lattice, or larger lattices, which soon 
becomes computationally prohibitive. 

An alternative approach to the calculation of the coefficients of the chiral Lagrangian 
has been suggested in ref. |]lOl- The basic observation there is that the separation 
between effective and non-effective degrees of freedom, which occurs in the continuum 
QCD theory, must be mirrored by an equivalent distinction in theory regularized on the 
lattice. Specifically, one can consider an effective Lagrangian, defined on the lattice in 
terms of an effective pseudoscalar field and external sources, which is equivalent to the 
fundamental QCD theory regularized on the lattice, for any value of the lattice spacing 
or bare coupling constant. 

This effective lattice theory can be then introduced as an intermediate step in the 
calculation of the continuum chiral Lagrangian. In order to derive the effective lattice 
Lagrangian, one can assume a sufficiently large set of couplings and fix the correspond- 
ing coefficients through a matching with an overcomplete set of expectation values. 
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computed both in the effective and the fundamental lattice theory. The several inter- 
actions, allowed in the lattice effective Lagrangian, can be organized as an expansion in 
terms of the distance of couplings, rather than in powers of the external momenta as in 
the continuum effective theory. Therefore, the lattice effective theory is not specifically 
a low-energy effective theory. This means that, in performing the matching between 
the two lattice theories, the existence of an infrared cut-off in the numerical simulation 
should not represent a problem any longer. 

A second advantage of this procedure is that the matching is performed between 
theories which are defined in the same regularization scheme, the four-dimensional 
lattice grid. For this reason, the finite ultra-violet cut-off effects, which affect the 
results of the numerical simulation performed with the fundamental lattice theory, can 
be kept better under control, because these effects are in principle exactly predicted 
and reproduced by the corresponding lattice effective theory. 

An important observation of ref. is that the lattice effective theory explicitly 
contains the collective fields which are responsible for the long distance behavior of 
the corresponding fundamental theory. For this reason, only short distance couplings, 
typically on the scale length of the order of the inverse rho mass, are expected to play 
a significant role in the effective theory. This observation then provides a criterion to 
select a finite number of possible interactions to be considered in the lattice effective 
theory. For the same reason, one expects that the determination of the corresponding 
coefficients should be feasible on a lattice of moderate size, thus allowing to achieve a 
better numerical accuracy in the calculation. 

Once the matching has been performed, and the lattice effective theory has been 
derived, one can consider the infinite volume limit of this theory and expand in increas- 
ing powers of the external momenta. The result of this expansion is the QCD chiral 
Lagrangian. In this way, the chiral coefficients can be calculated. 

A useful feature in this approach is that, in the limit of strong coupling and large 
number of colors A^, the integration @ of the non effective degrees of freedom in the 



QCD Lagrangian can be performed analytically on the lattice [|Tl| , |T2|. In this way, at 
the leading order of the strong coupling and large A^ expansion, the lattice effective 
theory can be exactly computed. 

This procedure has been followed by the authors of ref. in order to compute the 
coefficients of the QCD chiral Lagrangian in the strong coupling and large A^ limit. They 



4 



have considered a body centered hypercubical lattice, whose greater symmetry imphes 
that invariance under lattice transformations carries over to Euclidean invariance up 
to the terms of 0{p'^) in the continuum. Despite the strong coupling and large 
approximation, and an additional simplifying assumption of complete decoupling of 
the mesonic resonances, the final results for the chiral coefficients are in remarkable 
qualitative agreement with the experimental values. These results strongly encourage 
the attempt of a numerical calculation performed in the region of intermediate couplings, 
which is relevant for continuum QCD calculations. 

The feasibility of such a calculation is further investigated in this paper. With 



respect to ref. ||TOl, we consider the calculation on a standard hypercubical lattice, whose 
symmetry properties are well understood and four- dimensional Euclidean symmetry 
is known to be properly recovered in the continuum limit. Moreover, in this paper 



we study a fermionic action with the Wilson term |]13[, which serves to prevent the 
appearance of the undesired doubler fermionic species. A well known effect of the 
Wilson term on the lattice is the introduction of an additional source of chiral symmetry 
breaking, which persists even in the limit of vanishing bare quark masses. Main goal 
of this paper is a study of the effects of such a symmetry breaking from the specific 
point of view of the calculation of the coefficients of chiral Lagrangian, along the lines 
discussed above. 

In order to derive the low-energy effective theory for QCD, it is convenient to in- 
troduce a given set of external source in the fundamental Lagrangian In this way, 
the original global chiral symmetry of QCD can be formally promoted to an exact local 
invariance. On the lattice, because of the presence of the Wilson term in the fermionic 
action, it is necessary to introduce an additional set of external sources, which have no 
direct continuum correspondence. This step is discussed in sec. |], where we define all 
the external sources added to the action and the corresponding lattice operators which 
these sources are coupled with. 

In sec. H, we derive the corresponding lattice effective Lagrangian, at the leading 
order in the strong coupling and large expansion, by following the procedure of refs. 



10|-|12]. The main purpose of this calculation is to investigate the general structure of 



the lattice effective theory, in the presence of the Wilson term and the external sources. 
This general structure is a pure consequence of chiral invariance and discrete space- 
time symmetries, and must persist even in the region of intermediate lattice couplings. 
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Therefore, the lattice effective Lagrangian, as derived in the strong couphng hmit, can be 
assumed as a guidehne to define the corresponding theory in the region of intermediate 
couphngs. 

A primary effect of the presence of the Wilson term in the fermionic action is the 
appearance, in the lattice effective theory, of additional terms, which do not have any 
direct correspondence in the continuum chiral Lagrangian. The role of these terms is 
mainly to reproduce the effects of chiral symmetry breaking induced by the Wilson term. 
These effects are first investigated in sec. ^, where some typical physical quantities, like 
the pseudoscalar mass, the decay constant and the kaon Sx-parameter, are computed 
directly from the strong coupling lattice effective Lagrangian. 

The effects of chiral symmetry breaking are then further investigated in sec. ^, 
where we study the axial and vector Ward identities in the lattice theory. For the 
weak coupling region, these identities have been studied in ref. [Q. In that paper 
it was shown that, once the chiral limit of the theory is correctly identified and the 
lattice operators are properly renormalized, the predictions of continuum PCAC and 
current algebra are reproduced on the lattice up to vanishing cut-off effects, unaffected 
by the chiral breaking introduced by regularization. In this paper, we show that the 
same conclusion also applies to the lattice theory at the leading order of the strong 
coupling and large expansion. In this limit, the effects of the Wilson terms in the 
effective Lagrangian can be taken into account through a proper renormalization of the 
quark mass and the lattice operators. The correlation functions of these renormalized 
operators then satisfy, even in the strong coupling limit, the Ward identities of the 
continuum theory, as predicted by recovered chiral symmetry. 

In the effective continuum limit, that is in the limit of vanishing lattice bare cou- 
pling constant, the effects of chiral symmetry breaking induced by the Wilson term are 
expected to become completely negligible. However, it is well known that for typical 
values of couplings considered in current numerical simulations these effects are still 
quite relevant. For this reason, the role of the "Wilson terms", in the lattice effec- 
tive Lagrangian, cannot be neglected. On the other hand, the presence of these terms 
introduces additional difficulties. They significantly increase the number of couplings 
in the effective theory, they do not have a direct continuum limit in the QCD chiral 
Lagrangian and finally, because of their presence, the external sources in the lattice 
theory do not reduce straightforwardly to their direct continuum counterparts. 
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A procedure which allows to completely discard the Wilson terms in the lattice 
effective theory is then described in sec. ^. In this way, the general strategy discussed 
in this paper to derive the QCD chiral Lagrangian should be significantly simplified. 
The basic idea is to consider a new "renormalized" effective lattice Lagrangian, which 
reproduce the correlation functions of properly renormalized lattice operators in the 
fundamental theory. Because these correlation functions satisfy all the Ward identities 
predicted by continuum current algebra, the renormalized effective Lagrangian does not 
contain the Wilson terms at all. The way in which this mechanism works is illustrated 
through the discussion of few significative examples in the framework of the strong 
coupling effective theory, which exhibits the same chiral structure expected in the weak 
coupling limit. 

Finally, in sec. ^ we briefly discuss some aspects of the numerical calculation to 
be performed in the region of intermediate couplings of lattice QCD and present our 
conclusions. 



2 The chiral invariant Wilson action with external 
sources 

According to Gasser and Leutwyler a convenient preliminary step in the derivation 
of the QCD chiral Lagrangian consists in adding to the original QCD action a proper 
set of external sources. These sources are introduced in such a way that the resulting 
action becomes invariant with respect to combined local chiral transformations of the 
external sources and fundamental fields. In this section, we perform this step for QCD 
regularized on the lattice. But for the introduction of the Wilson term in the lattice 



action, our definitions of the external sources closely follow those of ref. ||10 



We consider the Wilson formulation of lattice QCD [|T^. The total action is: 



S=Sg + S^ (3) 

where Sg is the pure gauge action: 

^. = -^ETr(f/p + t/]>) (4) 
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Up being the plaquette variable and the sum over P running over all the oriented 
plaquettes on the lattice, and is the fermionic action: 



The quark fields, i/j and are n— component vectors in flavour space; are the lattice 
gauge variables, defined on the links, and m is the bare quark mass. For simplicity, we 
assume in this paper that the n flavours of quarks are degenerate in mass. 

In the limit of zero quark mass and vanishing Wilson term {m — r — 0), the lat- 
tice action is invariant with respect to global chiral transformations of the quark flelds: 
— > LiPl{x) and iPr{x) — > RiIjr{x), with L and R constant unitary matrices in 
flavour space. Formally, by adding to the fermionic action a proper set of external 
sources, this global invariance can be promoted to a local invariance, even in the pres- 
ence of quark masses and Wilson term. We consider a pair of "hnk-type" external 
sources, Wfi{x) and Z^{x), defined on the link between x and x + jl, and a pair of 
local sources, ^{x) and x(a;), defined on the site x. The sources W^{x) and Z^{x) are 
introduced to perform the parallel transport of the right- and left-handed quark fields 
respectively. The sources ^{x) and x(a;), which include the Wilson and mass term re- 
spectively, and are introduced to achieve chiral invariance of these terms as well. The 
lattice fermionic action, including the external sources, can be written as: 



(5) 



-|- (m + 4r) ^ \ijj{x)ijj{x) 



X 




+ ^{x + n) [rK{x + /x) + 7^] Ul,{x)Y^{x)iP{x) } + 



(6) 



-|- ^ ^ijj{x)J{x)iJj{x) 



X 



where the external sources, Yn{x), K{x) and J{x) are given by: 
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Wfi{x), Zf^{x), ^(x) and are matrices in the pure flavour space. 

It is easy to verify tfiat tfie total action, Sg + S^j, is now invariant with respect to 
local chiral transformations of the quark fields: 

^/'l(x) ^ V^l(x) , ipnix) R{x)ijR{x) (8) 

with L{x) and R{x) unitary matrices, provided the external sources are defined to 
transform according to: 



W^{x) ^ R{x)W^{x)R^{x + fl) 

r.) TMx 4- Ti,) 

(9) 



i{x) ^ R{x)i{x)V{x) 
X{x) R{x)xix) (x) 

The Wilson standard fermionic action is recovered from S^pj in eq. (|^) by con- 
sidering the limit: 

W^ix), Z,{x), ^ 1 ^^^^ 

In the following, we will refer to this limit as the limit of vanishing external sources. 
This limit will be also indicated with the symbol (. . .)q. 



By following ref. [|T0|, we express, in a non linear way, the external fields W^{x) and 
Z^{x) in terms of traceless hermitian vector and axial sources: 

= exp {-i + a^{x)]} ^^^^ 

Z^{x) = exp {-i [v^{x) - a^{x)]} 

where: 

v^{x) = v^ix) , a^{x) = a^ix) (12) 

and t'' are the generators of the SU (n) flavour group (A; = 1, . . . , n^ — 1). In the fermionic 
action @, the sources and are directly coupled to non-singlet vector and axial 
lattice currents. In the limit of vanishing external sources, the functional derivatives of 
the action with respect to and are given by: 

^ '''^ V'ix) , (-^] =Alix) (13) 



x 



/ 
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where 



+ iIj{x + ll)Ul{x) (7/, + r) t'' iIj{x) 



(14) 



and 



1/^*^(2;) is the extended lattice vector current, conserved for degenerate quark masses, 
and A^(a;) is the corresponding (non conserved) axial current. 

Similarly, we decompose the external fields ^{x) and x(x) in terms of hermitian 
scalar and pseudoscalar sources: 



^(x) = (7(x) + mix) 



X{x) = s[x) + iip{x) 



(16) 



Derivatives of the action with respect to these sources are given by: 



' 5S \ 

5a°'(x) 



S'^ix) 
-4rE"(x) 



' iSS \ 

Mix) Jo 
' iSS \ 



P'^ix) 



-ArWix) 



\"" V""/ / 

where S{x) and P{x) are the local scalar and pseudoscalar densities: 



s^ix) =i){x)ri){x) 



P"(x) = i){x)r'-i^i){x) 



(17) 



;i8) 



(a = 0, 1, . . . , — 1 and we use the notation — 1). The operators E(x) and n(x) are 
defined as: 



^"(^) = \Y. [i}{x)u^{x)t^i}{x + /i) + ij{x)ul{x - fi)rip{x - fi) 

n"(2^) = J E [i^{xWx)t-^,ij{x + /2) + V^(x)C/t(x - /x)r75V'(x - /2) 



(19) 



and are extended bilinears that, in the continuum limit, reduce to S{x) and P{x) 
respectively. 

Beside providing the formal invariance of the action with respect to local chiral 
transformations, the introduction of the external sources also offers a powerful tool 



10 



to compute the correlation functions of the lattice operators. Indeed, the partition 
function of the theory: 

Z{j) = f (dUdij(0) exp{-5(f/,^,^;j)} (20) 



is an explicit functional of the external sources (j = f^, a^, s,p, a, vr). We can then use 
eqs. (|T3|) and (|17D to compute, with standard techniques, the correlation functions of 
V^, A^, S and P as functional derivatives of Z{i) with respect to the proper external 
source^. Eventually, this approach becomes of practical utility once the partition 
function is expressed in terms of an effective Lagrangian. In the following, we discuss 
how this effective Lagrangian can, in fact, be calculated. 

3 The effective lattice Lagrangian in the strong cou- 
pling and large N limit 

In the limit of strong coupling and large number of colors A^, the high-energy degrees 



of freedom of lattice QCD can be analytically integrated out |Tl| , 0. After this inte- 
gration, the partition function turns out to be expressed in terms of an effective action, 
containing only mesonic fields and external sources. The purpose of this paper is to 
investigate the feasibility of performing the same integration, numerically, in the region 
of intermediate couplings relevant for simulations of lattice QCD. Still, the study of the 
strong coupling limit is an important preliminary step. The effective lattice Lagrangian, 
obtained in this limit, has a general structure, in terms of effective fields and external 
sources, that only follows from chiral invariance, charge conjugation and discrete lat- 
tice space-time symmetries. Consequently, the same structure must be preserved in 
the region of intermediate gauge coupling. In this region, one can consider an effective 
lattice Lagrangian of this form and determine its free parameters through a numerical 
matching with the full QCD lattice theory. 

The procedure to integrate out the high-energy degrees of freedom of lattice QCD, in 



the strong coupling and large limit, has been developed in refs. 12|. In ref. |10| the 



method has been implemented in the formalism of modern chiral Lagrangian introduced 



^Contrary to what happens in the continuum formahsm, however, here one has to be aware of the 
fact that the action S^j is not simply hnear in the external sources. Thus, the derivatives of the 
partition function might differ, from the corresponding correlation functions, for the contributions of 
additional contact terms. 
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by Gasser and Leutwyler 0. In this section, we apply the results of these papers to 
obtain the effective Lagrangian for lattice QCD with Wilson fermions and the set of 
external sources introduced in the previous section. 

At the leading order in the strong coupling and large expansion, the generating 
functional of eq. (120) can be expressed in terms of an effective action O, O: 



ZU) = I (dM) exp{-^e// iM;j)} (21) 
where Seff{M.;j) is given by: 

Seff{M;j) = -N |5:Tr[F(A,(x))]-ETr[lnA^(x)]+5:Tr[J(x)A<(a:)]| (22) 



and A4{x) is a matrix in the spin and flavour space representing an effective bosonic 
field. The function -F(A) in eq. (^) is defined as: 



1 + 



F{X) = 1 - VT^ + In [—^ j (23) 

and the matrix A^(x) is a function of the effective field and the external sources: 

A^(x) = -M{x) [rK{x) - 7^] Y^{x)M{x + Jl) [rK{x + fi) + YUx) (24) 



The external sources J{x), K{x) and Yf^{x) appearing in eqs. (p2D and (^J) have been 
defined in eqs. (|) and (|^). 

The functional integration in eq. (|2T|) is over a generalized contour in the complex 
matrix space, which can be parameterized by writing Ai in polar form: Ai = TZU, 
with TZ and U hermitian and unitary matrices respectively. The functional measure in 
eq. (pT]) is {dA4) = (dU), where (dU) is the Haar measure on the unitary group of U. 
The functional integral Z{j) is independent on the choice of TZ. 



In order to evaluate the partition function in perturbation theory, following refs. |TT 



r^l we look for a translationally invariant stationary point of the action Seff{Ai), taken 
in the the limit of vanishing external sources. The expansion of Ai around its vacuum 
expectation value will then generate the standard perturbative expansion and, for large 
N, this also corresponds to evaluate the partition function Z{j) by the saddle-point 
method. 
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The vacuum expectation value of Ai is found by choosing: 



Mo{x) = uol (25) 

and considering the effective action in the hmit of vanishing external sources. In this 
limit, the matrix A^(x) of eq. ( ^4|) reduces to: 

Xo = {l-r^)ul (26) 

The stationary point of S{A4) is then given by uq satisfying the equation: 

^(^-''^^'^o ^ (m + 4r) uo = l (27) 



whose positive solution is: 



-3 (m + 4r) + 4 V(m + 4r) + 7(1 -r^) 

Unim.r) = 4 (28) 

^ ' {m + 4rf + 16(l-r2) ^ ' 

For r = 1, one simply has: 

Uo{m,l) = ^— (29) 
m + 4 

The matrix Ai can be now expressed in terms of fields with vanishing vacuum 
expectation value. A convenient general form is: 

M{x) = uo exp[iS{x) +iV{x)-f5 + iV^,{x)-f^ + iA^,{x)'y^'y5 + iT^^{x)a^u] (30) 

where S, V, V^, and 7^,^ are matrices in the pure flavour space. In the lattice 
effective theory these fields, which carry integer spin, represent the mesonic excitations 
with the corresponding quantum numbers. 



The lattice effective action of eq. (22) differs from the standard continuum chiral 
Lagrangian in that the latter is only expressed in terms of external sources and pseu- 
doscalar fields. Thus, in order to derive an effective lattice action corresponding to this 
continuum Lagrangian, the contributions of the scalar, vector, axial and tensor com- 



ponents of in eq. (|30D , should be integrated out. In principle, this could be done 
explicitly. However, in the context of the strong coupling approximation considered so 
far we are more interested in deriving the general structure of the effective Lagrangian 
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than in computing its specific expression. Therefore, also in the attempt to keep dis- 
cussion simple, we will not perform such an integration. Rather, we assume a complete 
decoupling of the higher resonances (corresponding to the infinite mass limit) by simply 
neglecting their contributions. That is, we constrain j\4{x) to the form: 



M{x) = uq exp [iV{x) 75] = uo 



(31) 



where U (x) = exp [iP(a;)] is a unitary matrix in the flavour space. In the continuum 
limit, U{x) would correspond to the standard effective field entering in the QCD chiral 
Lagrangian. 

An explicit expression for the effective action, in terms of the field U, is obtained 
by Taylor expanding the function -F(A) of eq. (p^) around the stationary point A = Aq. 
The resulting action has then the form of an infinite series: 



k=l 



where: 



and 



SmiU) 



Sk{U) = 
NC„ 



Y^TT[X,{x)-Xof 



^Tr M(x)?7(x) + x{x)U^{x] 



The coefficients Cm and Ck {k = 1,2, . . .) are given by: 

4u2'=F«(Ao) 



Cm = 4^0 



Ck 



k\ 



(32) 

(33) 
(34) 

(35) 



and are of order one in the large N expansion. For practical calculations, the spin- 
flavour matrix A^(x) — Aq can be conveniently written in the form: 



A^(x) - Ao = - ^ [- U{x) (A^f/(x))^ (1 + 75) - U\x) {A^U{x)) (1 - 75) + 
+ r^Uix) {np{x)y (1 + 75) + r2f/t(x) (n^Uix)) (1 - 75) + 
+ rUix) (S^t/(x)) ^ (1 + 75) 7m + rWi^) {^^Uix)) (1 - 75) 1, 



(36) 
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where A^U is: 

A^t/(x) = W^{x)U{x + V)Zl{x) - U{x) (37) 
and the quantities ^^'^U and T^^'^U are defined as: 

n^U{x) = W^{x)^{x + i2)W{x + Jl)Wl{x)^{x) - U{x) 

n^^Uix) = ax)Z,ix)Wix + + f^)Zlix) - U{x) ^^g^ 
^p{x) = Z^{x)W{x + n)W^^{x)ax) - Z^{x)^^{x + fi)U{x + fi)Zl{x) 
^?U{x) = ^{x)Z^{x)W{x + n)Wl{x) - W^{x)U{x + fi)^^{x + fi)Wl{x) 

Equation (|3^) is the effective chiral Lagrangian for lattice QCD with Wilson fermions 
at the leading order in the strong coupling and large expansion. In the next section 
we will study the predictions of this Lagrangian in some more detail. Here we conclude 
with some important remarks. 

i) As the original fermionic theory, the effective Lagrangian is invariant with respect 
to local chiral transformations. This is easily shown to be the case when we assume 
that the field U (x) transforms according to: 

U{x) ^ R{x)U{x)L\x) (39) 

and the external sources W^, Z^, C, and x have the transformation properties defined 
in eq. (^. Notice that the quantity A^U{x), introduced in eq. (|37|) , transforms as the 
matrix U{x) and, in the limit of vanishing external sources, it reduces to the standard 
discrete derivative of U: 

(A^[/(x))o = V^f/(x) = U{x + U{x) (40) 

Thus, A^f/ represents a possible definition of the lattice chiral covariant derivative of 
U and, in the continuum limit, it reduces in fact to the continuum covariant derivative 

ii) In the limit of vanishing Wilson term (r = 0), the quantity (A^ — Aq) in eq. ( p6D 
is proportional to the covariant derivative A^U{x). Therefore, it vanishes for vanishing 
external sources and for constant field U{x). In this limit, the series expansion of 
the effective action in eq. (|32|) corresponds to an expansion in powers of the external 
momenta p, each term Sk being at least of order p^. As in the continuum case, this 
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feature is a consequence of chiral symmetry, and is what makes the continuum QCD 
chiral Lagrangian a useful tool for the study of low-energy strong interactions. However, 
due to the chiral breaking of the Wilson term, this feature is lost on the lattice. In 
particular, contributions of the order of enter in each of the infinite terms Sk of the 
effective action. Thus in this case, even in the low-energy limit, each of these terms 
should be in principle considered to perform calculations. Nevertheless, as we will show 
in the next section, in the large N limit one can use the standard perturbative technique, 
so that the lattice effective action of eq. (|3^ ) is still useful in practice. 

iii) The effective lattice action of eq. ( ^2]) is not specifically a low-energy effective 
theory. Indeed, but for our deliberate neglecting of the contributions of the higher reso- 
nances, no low-energy approximation has been introduced to derive it. Rather, eq. ( ^2|) 
defines an effective theory describing, in the strong coupling and large limit, the 
strong interactions among pseudoscalar mesons at arbitrary energy. In the interme- 
diate coupling region, a corresponding effective Lagrangian, with the same symmetry 
properties, can also be considered. The free parameters entering in this Lagrangian, 
including those coming from integrating out the higher resonances, can be eventually 
computed by matching with a numerical simulation within the full lattice theory. In 
this approach, the infra-red cut-off introduced by the finite lattice size in numerical 
simulations does not represent a problem any longer, since no low-energy expansion has 
been performed in deriving the effective theory. 

iv) The most evident effect of the Wilson term in the lattice effective action is the 
appearance of additional couplings. For r = 0, the effective action of eq. ( ^21) reduces 
to the form: 



Sd{U) = -N ETrU^i 



2C, 



(A,f/)t (A,f/) 



(41) 



which, but for the obvious effects of discretization and lattice breaking of Lorentz in- 
variance, has a direct correspondence with the continuum chiral Lagrangian. 

The presence of the Wilson term in the original lattice action gives rise to additional 
interactions in the effective theory, of the kind: 



Sw{U)=Nr^ ETr -Ci 



+ 



+ C2 
-u 



(42) 
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These terms are explicit sources of chiral symmetry breaking, since in the hmit of van- 
ishing Wilson sources, ^ ^ 1, they reduce to interactions of the form f/^, (A^U) (A^U), 
. . . . In the continuum limit {g 0), the coefficients of these terms are expected to 
vanish. However, in the region of intermediate couplings relevant for current lattice 
simulations, their effects are still not negligible. These terms significantly increase the 
number of allowed couplings in the lattice effective theory and, furthermore, they do not 
have a direct correspondence in the continuum QCD chiral Lagrangian. In the following 
of this paper we will discuss the way to overcome this problem. It will be shown that 
an effective lattice theory, without the Wilson terms, can be in fact considered, and this 
theory reproduces the values of properly renormalized lattice correlation functions. 

4 Correlation functions and matrix elements from 
the effective Lagrangian 

In this section we use the effective lattice Lagrangian, derived in eq. (p^), to calculate 
explicitly some simple physical quantities. From a phenomenological point of view, 
the results are significantly affected by the approximations introduced so far: strong 
coupling limit, large N expansion and assumption of complete decoupling of the higher 
resonance. On the other hand, these calculations achieve a twofold goal: first, they offer 
us the opportunity to illustrate the perturbative technique, as applied to the lattice 
effective Lagrangian. Possibly, the same technique can be used for determining the 
coefficients of the effective Lagrangian in the intermediate coupling region. In addition, 
the analysis which follows, combined with the study of chiral Ward identities in the 
next section, will point out some important effects of the several Wilson terms entering 
in the lattice effective theory. 

As discussed in sec. ^ the correlation functions of the lattice operators can be 
calculated by differentiating the partition function Z{j) with respect to the external 
sources. The first convenient step in order to do that, consists in expressing the effective 
field U{x) in terms of — 1 independent real fields, representing the true degrees of 
freedom of the effective theory. We choose the following parameterization: 

U{x) = exp {2 i 0(x)/Fo} (43) 

where 4>{x) = 4'k{x)tk and Fq is a coefficient introduced to conventionally normalize 
the pseudoscalar field. 
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Since the derivatives of the partition functions must be eventually computed in the 
limit of vanishing external sources, it is also convenient to perform, from the beginning, 
the following change of variables for the external scalar sources: 

soix) = {m + Ar) + s'q{x) , aoix) = 1 + a'^ix) (44) 

In this way, the new sources Sq(x) and c'^lx) (we will omit the primes in the following) 
do actually vanish in the limit of vanishing external sources. 

The next step consists in expanding the effective action (p2D in increasing powers 
of the pseudoscalar field and the external sources. Since Aq represents the vacuum 
expectation value of the matrix A^(x), the difference A^(x) — Aq is, by construction, of 
the order of the field and the external sources. Thus, although the action (^) has 
the form of an infinite series in powers of (A^(x) — Aq), the above expansion can be 
systematically performed. For convenience, we fix the parameter Fq of eq. (H3|) in such 



a way that the the kinetic term of the effective Lagrangian for the real scalar field has 
the standard normalization, i.e. (1/2) (V^0)^. In this way we obtain: 

= 2N (l + r^) \Ci + 2C2 (l - r^)] (45) 



where Ci and C2 are the coefficients defined in eq. (pSj) . The effective action 5*0, as a 



function of the field 4> and the external sources, can be written in the form: 

^<^=^r+^?^+^r+... (46) 

where the generic term S*^^ contains the external sources at the power k. The action 
S^^ represents therefore the effective lattice action in the limit of vanishing external 
sources. It is given by: 



(47) 

+ ^43 K ( V^03) (y^^) + ^44 + . . . } 



For simplicity, all the fiavour indices in the above expression have been omitted. Thus, 
for instance, (V^0)^ stands for (V^0i) (V^0i) and K(j)'^ for Kijki(f)i(j)j(j)k4>i- The tensor 
Kijki is given by: 

Kijki = —SijSki + -dijsdkis (48) 
n 2 
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and the structure constants dijk, as well as the constants fijk to be introduced below, 
are defined by the algebra: 

1 . 1 

The coefficient of the 0^ term in eq. (^) (we are anticipating that this coefficient 
represents in fact the pseudoscalar meson mass) is given by: 

^' = -pf^(^ + 4r)-(r^ (50) 
and the couplings dij of the several 0^ terms are: 

(1 + 6r2 + r^) ^2 , _ (1 - ^1 - 64r2A2 

W41 — o 1 U42 — o 

2N{l + r^fAl 8N{l + r^fAl 

[1 + r^) Ai - 24r2 (3 + r^) A2 . _ 8r^ {Ai - Ur'^Ai] 

dd 



(51) 



^43 — TTTTV-, ox2 .0 ' "44 



where: 



6Ar(l + r2)Mf 3Ar(l + r2)M2 



Ai = Ci + 2C2(l-r2) 

A2 = C2 + 3^3(1 -r2) + 2^4(1 -r^)' 



(52) 



Since, at the leading order, the effective action S^ffiU) is proportional to the number 
of colors A^, the perturbative expansion in terms of the scalar field (f) corresponds in 
fact to an expansion in increasing powers of 1/A^. The couplings of the cfP' terms are 
by construction of order 1 and the couplings of the terms are of order 1/N^~'^ 
(see eq. (|5TD). Thus, in the large N limit, the lattice effective Lagrangian can be 
treated in perturbation theory. Furthermore, it is important to remember that, in this 
context, only the lowest (non trivial) perturbative order should be considered in each 
particular calculation, since contributions suppressed by additional powers of 1/A^ have 
been already neglected, from the very beginning, in the calculation of the effective 
Lagrangian itself. 

At the leading order in perturbation theory, the propagator of the scalar field 
satisfies the equation: 

(V^-M2)A(x,y) = -5.,, (53) 
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where represents the (symmetric) discretized version of the four-dimensional Lapla- 
cian operator: 

V'fix) ^ Yl Ifi^ + 2/(x) + fix - n)] (54) 

The solution of eq. (^) is the standard lattice propagator of the scalar field: 

1 p-ip{x-y) 
A(x, 2/) = 77 y — (55) 

where V = L^LyLzLt is the lattice volume and the sum is extended over all the lattice 
momenta, pi = 27ik/Li. For later use we also give the expression of the propagator as 
a function of the time t and the spatial momentum p. In the infinite volume limit, and 
for t > 0, one finds: 

Ait,p) = y A(t,f;0) e-^^^-'=-^^ (56) 
^ 2 smh 



where Et, is given by: 

,1/2 

^ / rt; \ , / /I/; _ 1 

= 2arcsinh 



■ . 2 fPi\ ^ . ^2(M.. 

sm — + smh ' 



2 



(57) 



M^ = 2 arcsinh (^-^ ) (58) 



In the continuum limit, — > M^, and E^, —>■ J pf + M^. 

/-I \ * 

The next term, , of the lattice effective action in eq. (|46| ) is the one linear in the 
external sources. This term can be conveniently written in the form: 

= V |so(x)S'o(x) + (To (a;) So (a;) +pk{x)Pk{x) + 7[k{x)Ilk{x) + 

X ^ . (59) 

+a^^{x)A^^{x)+vf,{x)V;:{x) + ...} 

where the "effective" operators Oi are given by: 



^o( 


x) 


— ds + dsi 


i)l{x) + ... 


So 


» 




<j)l{x) + ... 


Pk( 


x) 


= dp(j)kix) 


+ ... 


flk 


[x) 


= dj,i(j)k{x) 


+ d^2V^(f)k{x) + 


Ik 


(x) 


= daVf^(j)k{ 


x) + ... 


yk 


(x) 


= d^ fijk (pi 


(x)(V^0,(x)) + .. 



(60) 
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with coefficients: 



dp = -NCjFo , da = -Fo/il + r') , d, = 1 

(61) 

d^ = SnNCir'^ , d^i = -8rV (1 + r^) 

d^i = 8r2Fo/(l + r2) , d^2 = {1 + r^) 

In the effective theory, the operators Oj can be used to calculate directly the correlation 
functions of the corresponding operators in the original theory. 

When one needs to compute the correlation functions of higher powers of the lattice 
operators (like for instance AjJ then additional contributions might come from terms 
containing higher powers of the external sources (like and a^0^). Examples of such 
terms, that will be considered in the study of the lattice Ward identities, are found in 
the next term of the effective action (HSI): 



X 

+ d^^ fijk [(V^0i(x)) 7ij{x) - 0i(x) (V^7rj(x))] v^{x) + 
+ dav fijk a'^{x)vj^{x) [0fc(x + /i) + 0fc(a;)] + 

~l~ dyy fijkflmk 

vl^{x)v^^{x)(j)j{x)(j)^{x + fl) + ...} 



(62) 



where: 



rf^^i = 2NC,r^ - 2r'F^/ (1 + r^) , d^^2 = r^F^/8 (1 + r') 

daa = -da. = F^/2{l+r') , d^, = Fo/2 {1 + v') (63) 

duy 

In the effective theory, the correlation functions are computed by expressing the 
partition function Z{j) in terms of a functional integral on the fields (f). Up to higher 
terms in 1/A^, the functional measure simply reduces to (c/0): 

ZU) = J m exp{-540;j)} (64) 

All the correlation functions of interest can be then obtained by subsequent derivatives 
of Z{i) with respect to the external sources. We conclude this section by discussing a 
few interesting examples. 
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The quark condensate 

The simplest relevant case is the calculation of the quark condensate {ipip). This quan- 
tity vanishes to all orders in standard perturbation theory and a value of {ipip) different 
from zero is expected to be the signal of spontaneous chiral symmetry breaking. In the 
effective theory, at the leading order of the strong coupling and large expansion, the 
absolute value of the condensate for a single flavour of quarks is given by: 

Thus Mo, the vacuum expectation value of the effective field A^(x), basically represents 
the quark condensate. 

The pion mass and the critical quark mass 

In the effective theory, the correlation function of two pseudoscalar densities is given 
by: 

where A(a:, y) is the propagator of the 0-field of eq. ([551). At this order, the pseudoscalar 
correlation function only receives contribution from the propagation of a single pseu- 
doscalar meson. Thus, A(x, y) is just the pion propagator and the pole of this function, 
M^, defines the pion mass. From eq. (|50|) we can write this mass in the form: 

2NC„ 



= — — 2^ (m — m(m, r)) (67) 



where m{m, r) is the function: 

^^"^'^)=(rT^^)A^-'^ (68) 

The vanishing of the pion mass identifies the chiral limit of the theory. The corre- 
sponding critical value of the quark mass, m-c, is therefore defined by the equation: 

TTLc = m{mc, r) (69) 

In the limit r = 0, the function m identically vanishes and the critical quark mass is 
rrir = 0. For finite values of r, the Wilson term induces an additive renormalization of 
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the quark mass. In particular, for r = 1 one finds from eqs. (^) and (^) that the pion 
mass vanishes for rric = —2, corresponding to the critical value kc = 1/4: of the Wilson 
hopping parameter. This is indeed a well known result of strong coupling lattice QCD0. 

In numerical simulations, it is often convenient to compute the correlation functions 
as a function of the time t and the spatial momentum p. In the effective theory (in 
the infinite volume limit) these functions can be calculated by using eq. (|56|) . Thus, for 
instance, for the pseudoscalar correlation function we find: 

Cip(^,Pl = T.{P'i^)P'^me-''' = (70) 
- 2 smh 

with given by eq. (0). It is interesting to notice that the energy dependence of the 
correlation function predicted by the effective Lagrangian differs, by lattice discretiza- 
tion effects, from the corresponding continuum form: 



where E^, = (p^ + M^)^^^. A behaviour as given in eq. (ffOD has been observed in 



numerical QCD simulations [|I^]. Therefore, this is an example of pure lattice artifacts, 
that are actually found in numerical calculations, and whose existence is predicted by 
the lattice effective Lagrangian, even in the strong coupling approximation considered 
so far. It is likely that, in the numerical determination of the coefficients of the effective 
Lagrangian in the intermediate coupling region, the possibility of taking into account 
such lattice artifacts will help in better controlling the associated systematic errors. 
This is another advantage of introducing the effective Lagrangian on the lattice as an 
intermediate step in the calculations. 

The pseudoscalar decay constant 

The decay constant of pseudoscalar mesons can be computed from the following corre- 
lation function: 

<P(^,P) = T.{^ii^)P''me-''^-' = -6^' (l - e-^A e-^^* (72) 

^ 2smh£/7r ^ / 



^It can be proved, in fact, that the result at r = 1 is true for any N , and applies equally well to 
iV = 3 111. 
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In general, two-point correlation functions in momentum space can be used to extract 
the values of the matrix elements of lattice operators between the vacuum and a single 
particle state. For instance, the correlation function of any two operators Oa and O^, 
carrying the quantum number of the pion field, are expected to have in the large time 
limit the following behaviour^ 

2 sinh 

At the lowest order in the strong coupling effective theory the large time limit is not 
a necessary requirement, since the correlation functions are always of the form ([75|). 
Thus, from eqs. (|7Dp and (^), we can obtain (up to irrelevant complex phases) the 



values of the following matrix elements: 

(0|P^(0)|7r'=) = 5^X ' (0|4(0)|7r'^) = -5^X(l-e-^-) (76) 

The last equahty provides the value of the pseudoscalar decay constant F^, in the strong 
coupling limit. By considering the continuum definition (0|Ao(0)|7r*'') = S^^Ft^Et^, we 
find that, up to lattice artifacts: 

F^ = -da = Fo/{l + r^) (77) 

From eqs. (0), (0) and ([77|), we can also derive a relation between and the quark 
condensate: 

(l + r^y Fl Ml = 2{m-m) {^ijj) (78) 

In the next section, the same relation will be derived directly from the axial Ward 
identity of the lattice theory. In that context the meaning of the factors (1 + r^), 
entering in eqs. ( [77| ) and ( |78D because of the Wilson term, will be explicitly clarified. 

■^For convenience we choose here, for the one-particle states, the normahzation: 

{Mip)\M{p')) = 2 sinh Em6pp> (73) 
instead of the standard covariant normahzation: 

{M{p)\M{p')) ^ 2EmSpp' (74) 
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The form factor of semileptonic pion decay 

In the hmit of exact isotopic invariance, the semileptonic decays of charged pions, 
TT^ TT^e~^i>e, are described in terms of a single invariant form factor, This form 

factor parameterizes the matrix element of the weak vector current between external 
pion states: 

(7r°(p2)|«7M^|vr+(pi)) = V2U{q') {pi+p^}, (79) 

<i^i = {pi — P2)fj. is the momentum transferred in the decay and the form factor is 
normalized at = by current conservation: /+(0) = 1. 

The matrix element of eq. can be computed by considering the following three- 
point correlation function: 

g-E2ty + Elt:^ 

4 sinh El sinh E2 

where the last equality holds for tx < < ty. An explicit calculation of this function 
in the effective theory shows that the matrix element of the weak current between two 
external pion states has the form: 

{^'iP2)K\0)\7r\pi)) = -id^f.ui {e'^'^ - e-^^^") (81) 

Thus, up to lattice artifacts, we find that the vector form factor for pion semileptonic 
decays is just given by: 

/+(g^)=4 = l (82) 

This result has likely a simple explanation. In reality, the form factor f+{q^) turns 
out to be well described by a phenomenological vector meson dominance model, where 
it is expected to have the form: f+{q^) = rrip/ {nip — g^), with uip the rho meson mass. 
In this case, in the limit of complete decoupling of the rho meson as the one considered 
so far, the form factor becomes completely independent on and reduces to 1 for 
current conservation. In a typical lattice simulation, the value of /rn^p is roughly of 
the order of 1. Thus, in that context, the assumption of complete decoupling of the rho 
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meson is certainly unreliable. What we expect to happen is that the effect of the rho 
and other resonances will manifest itself through the structure of the couplings in the 
effective lattice Lagrangian, which, rather than being limited to the nearest neighbors 
as in the strong coupling approximation, will fall off with a range ~ rrip^. It should be 
also noticed that the equality = 1, obtained for the coefficient of the effective theory, 
is a consequence of the lattice vector current conservation. In the next section we will 
find in fact that the same equality can be also obtained directly from the vector Ward 
identities of the lattice theory. 

The 5/^-parameter 

The calculation of the iJ/^-parameter in the effective theory shows that some chiral 



symmetry breaking effects, induced on the lattice by the Wilson term |T^, |T6[ and 
actually observed in numerical simulations (see e.g. ref. do not appear in the 

strong coupling and large N limit. This seems to be a quite general property of the 
theory in this limit: that is, the chiral symmetry is not broken to all the extents it 
could possibly be. The same feature will be illustrated by the study of the lattice Ward 
identities in the next section. 

The parameter parameterizes the matrix element (K^\0^^^'^\K^) relevant for 
K^—K^ mixing in weak interactions. O^^^^ is the four-fermion operator (s7^(l — 75)^)^ 
entering in the effective AS* = 2 weak Hamiltonian. For external states of kaons at rest, 
the parameter Bk is defined as: 

(K°|0^^=2|ir°) = ^-^BkFImI (83) 
This definition is such that in the vacuum saturation approximation: 

so that, in this limit, B^ = 1 for = 3. 

On the lattice, with Wilson fermions, the matrix element between external states of 
arbitrary momenta can be conveniently parameterized in the form: 

(Z°(p2)|0^^=2|K°(pi)) = a + /5Mi + 7 (Pi ■ P2) + . . . (85) 

where the dots stand for higher order terms in M^- and (pi ■ ^2)- The parameters a and 
j3 are lattice artifacts that should vanish in the true continuum limit, g ^ Q. The origin 



Bk = -Ai + \-A (84) 
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of these terms is the mixing of the lattice operator O^^ ^ with operators of different 



chirahty , and it is only due to the presence of the Wilson term in the lattice action. 

The calculation of the matrix element ( pSf ) in the effective theory proceeds through 
the evaluation of the corresponding three-point correlation function, by using the tech- 
nique discussed above. In the large N limit, one finds that the matrix element only 
receives contribution from the part of the (V" — A)"^ lattice operator. The matrix 
elements of the VA and AV operators vanish identically for parity invariance and the 
remaining contribution, V"^, is suppressed by an additional power of 1/N . Technically, 
this follows from the fact that the relevant coefficients in the effective Lagrangian are 
da ~ \fN and = 1 (see eq. (pl|)). In addition, we find that in the large N limit the 
matrix element is saturated by the vacuum insertion and it is given by: 

(Z°(p2)|0^^=2|ir°(pi)) =Adl (e'P^^ - l) (e-'P^^ - l) (86) 



Up to discretization effects, eq. (|86| ) corresponds to eq. (|^) with a = /3 = 0, as 
would be predicted by unbroken chiral symmetry. Therefore, at the leading order of 
strong coupling and large expansion, the mixing of O^^^"^ with operators of different 
chirality does not occur. We also find in this limit 7 = 4(i^ = 4F|;, corresponding to: 

Bk = \ (87) 

According to eq. (|8^) , this is the result predicted by the vacuum insertion approximation 
in the large N limit. 

Although the 1/A^ corrections to eq. ( p7D cannot be consistently computed within 
the approximations considered so far, it may be useful to observe that 1/A^ suppressed 
contributions to the parameter a in eq. ( ^5]) in fact exist. They come for instance from 
the Wilson term: 

- NCir^ Tr [a^.{x)(P{x + /i)a^(x)0(x) + (a^ v^)] (88) 
x,^l 

which can be obtained by expanding the effective action Si of eq. (|32D|^. Other 1/A^ 
suppressed contributions to a are also found at the next order in perturbation theory, 
coming from the pure 0^ term of the effective action (see eq. (0)). Since the coefficient 
(^44 of this term is proportional to r^, it represents again a pure effect of the Wilson 
chiral symmetry breaking. 

"'The same combination (^) is also found in the symmetric term (A^t/)^ (A^t/) of the effective 
action, but now with a relative minus sign between the vector and axial part. Thus, in this case, the 
contributions of the two terms cancel in the final result. 
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5 Chiral invar iance and Ward identities 



The chiral invariance of the Wilson action, in the presence of external sources, implies 
the existence of a full set of vector and axial Ward identities (WI). These identities 
have been studied in ref. |n] to show that, in the continuum limit, they reduce to the 
known WI's of QCD. Precisely, once the chiral limit of the theory is correctly identified 
and the lattice operators are properly renormalized, the WI's on the lattice reproduce 
the relations of continuum current algebra. 

In this section, we will first verify that the lattice WI's are indeed satisfied, by 
explicitly computing the relevant correlation functions in the effective theory. This is 
nothing more than a useful check of the calculations performed so far. Then, however, 
by using these results, we will show that the relations of continuum current algebra are 
exactly reproduced by the lattice theory even in the (leading order) strong coupling 
and large expansion considered so far. In a sense, this is quite a unexpected result. 
According to ref. the current algebra relations are expected to be recovered only 
in the continuum limit, when the contribution of higher dimensional operators to the 
correlation functions of the WI's can be exactly neglected. In contrast, the strong cou- 
pling approximation represents rather the opposite limit, in which such contributions, 
if any, should appear as genuine and finite corrections. 

For the purpose of this paper, the possibility of recovering exact chiral symmetry 
with the Wilson action in the strong coupling approximation is a very useful result. 
In fact we can define, even in the strong coupling limit, "renormalized" operators that 
satisfy all the WI's of the continuum theory. Then in the next section, by consider- 
ing some specific examples, we will be able to show, analytically, that the correlation 
functions of these renormalized operators are completely reproduced by a lattice ef- 
fective Lagrangian which does not contain the Wilson terms at all. Thus, each term 
of this Lagrangian, but for discretization effects, has a well defined correspondence in 
the continuum effective theory. The same result must be valid in the scaling region of 
couplings relevant for QCD simulations. In that case, this effective lattice Lagrangian 
can be eventually expanded in powers of the lattice spacing and external momenta to 
compute the coefficients of the continuum QCD chiral Lagrangian. 

In order to study the lattice WI's in the framework of the effective theory, it is 
convenient to derive them in the form of equalities among derivatives of the partition 
function with respect to the external sources (see footnote 1). 



28 



Let us consider the set of infinitesimal vector and axial transformations generated, 
from eqs. (|^) and in the limit in which: 



R(x) ~ 1 + i \av(x) + aA(x)] 

(89) 

L{x) ~ 1 + i [a;v'(x) — ^^(a;)] 

where av,A{x) = ay^(x)t'^ are infinitesimal quantities. In this limit, the chiral trans- 
formations of the quark fields in eq. (§) reduce to: 

ip{x) ip'{x) = ip{x) + i [avix) + aA(a;)75] ip{x) 

— — / — — y^^) 

ip{x) ^ {x) = ii{x) - iip{x) [av{x) - 750:^(3;)] 

Similarly, from eq. (j^), we can derive the transformations of the external sources. For 
the vector and axial sources they have the form: 

5Vf,{x) = V^av{x) + - [av{x + fi) + av{x),Vf,{x)] + 
i 



+ - [aA{x + f2) + aA{x), a^(x)] + . . . 
Sa^{x) = V^aA{x) + - [av{x + fi) + av{x), a^{x)] + 
+ - [aA{x + + aA{x),Vf,{x)] + ... 



(91) 



where the dots indicate terms of higher order in the external sources. For the scalar 
and pseudoscalar sources these transformations, by being local, have the same form 
expected in the continuum theory: 

6a (x) = i[av{x),(7{x)] — {aA{x),n{x)} 
5'k{x) = i[avix),-K{x)] + {aAix),a{x)} 

(92) 

Ss{x) = i[av{x), s{x)] — {aAix),p{x)} 
Sp{x) = i[avix),p{x)] + {aAix),s{x)} 

Let us then perform in the partition function Z{j) of the full lattice theory the change 

of variables defined by eq. (pH). With respect to this change of variables the measure 

in the functional integral is invariant^. On the other hand, by definition, a change of 

^In the more general case of U{n) (g) U{n) flavour transformations the invariance of the measure is 
lost and the consequent variation generates the quantum anomaly. 
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variables does not affect the value of the integral. Thus, by considering that the lattice 
action is invariant with respect to combined transformations of the quark fields and 
external sources, it is easy to proof that the partition function itself is invariant with 
respect to transformations of the external sources only. Indeed, by denoting with a 
prime the transformed quantities, one has: 

= f [dUdij(hp) exp{-S{u,ij,t,j)} = 

[dUdtfjchp) exp { -5 (u, tp', j') } = (93) 
(dUdij'dW) exp {-S [U, W; f) } = Z{f) 
Formally, we can express this invariance in infinitesimal form: 



Sjaix) = 



(94) 



where Sj{x) are the variations of the sources given in eqs. (|9TD and (|92D. The whole 
set of vector and axial WI's are finally obtained by further differentiating the above 
equation with respect to the external sources and taking the limit of vanishing external 
sources. For example, for the two- and three-point correlations functions, the generic 
WI's have the form: 



EE 



and 



EE 



^Sjaix) Sjb{y)) \Sjc{z) 







(95) 



+ 



,5ja{x) 5jc{z)) \5jbiy) J ' \5jaix)J \5jc{z) 5]b{y), 



+ 



(6m \ / 5^ja{x) 



(96) 







respectively. 

To begin the analysis of the WI's in the effective theory, let us consider the simple 
case of the two-point axial-pseudoscalar correlation function. The corresponding WI is 
obtained from eq. (^) by choosing i^{y) = p'{y)- The result reads: 



1 



6^Z 



Z 6a^^{x)5pi{y) 
' 1 



2 (m + 4r) 



1 



5''Z 



+ 2 



' 
5^Z 



Z 6T\^{x)6p{y) ) „ n 



- 5^^ 6^ 



Z 6p'^{x)6p^{y) ^ 
( 1 5Z 



(97) 



x,y 



Z5s^{x) 
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dp (4i + d^2M^) 6'^^ A{x, y) - d,dn2 S'^6,,y (99) 
dl5'^^{x,y) , (It^) =-ds 



where Vf; denotes the left discrete derivative: 

V^^f{x)^f{x)-f{x-n) (98) 

and, we recall, n in eq. (pTj) is the number of flavours. 

In the effective theory, the derivatives of the partition function with respect to the 
external sources can be easily calculated . For the terms entering in eq. (|97|) we find, 
at the lowest order: 

(^H(»)o ^ '"'^^^^'^^'^^^'^^ - '^'^'^''-^ 

\Z5Ti^{x)5jP{y))^ 
fl \ =d^5^^Mxy) 

Now, for the WI (P7| ) to be satisfied, the terms proportional to A{x,y) and those 
proportional to the 6x,y must be equal on both sides of the equation. From these 
requirements we find two relations among the several coefficients di: 

{da - 2^2) = 2 (m + 4r) dp + 2d^i (100) 

and 

{da - 2d^2) dp = -djn (101) 

By using the explicit values of the coefficients, (eqs. ( pT| ) and (pT])), it is easy to verify 
that both the relations ( |100| ) and ( |101| ) are indeed satisfied. Thus, in the effective 
theory, the axial WI (|97D is identically satisfied as well. 

To extend this check, we have considered a larger set of both vector and axial WI's, 
for two- and three-point correlation functions. In each case, one or more relations among 
the coefficients di are found, each of the kind of eq. ( |101|) . Examples of such relations, 
as derived from the axial WI's for the correlation functions (AH), {AA), {ASP) and 
{AAV), are: 

{da - 2(i^2) d^i = -d^/n + M^^i , {da - 2d^2) dn2 = -4:d^n2 

{da - 2d^2) da = 2daa , {da - 2d^2) 4i = dp (102) 

{da 2c?7r2) d^ 2dau 2d^y , {da 2c?^2) dy da 2d-j^y 
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Similarly, by considering the vector WI's for (V^V"), (VPP), (VPU) and (VAA), we 



find: 



dv 



dr., 2 6^7) 



d-K2 



da 



da/2 (103) 

All the relations in eqs. (|102| ) and (|103|) are satisfied regardless the specific values of 
the coefficients Cm, Ci, C2, ... entering in the effective Lagrangian. In other words, 
eqs. ( |102|) and (|103| ) are a pure consequence of chiral symmetry and the fact that these 
relations are satisfied is basically just a useful check of the calculations performed so 
far. 

It is useful to rewrite eq. ( |10CI| ) in the following form: 



where we have defined: 



{da - 2d^2) Ml = 2[m -ffi (m, r)] dp 



m (m, r) = — — 4r 
dp 



(104) 



(105) 



We easily recognize that the identity (|104|) is just eq. (^), that fixes the ratio between 
the square of the pion mass and the quark mass. In particular, the function m, defined 



in eq. (|105|) , is exactly the same quantity previously introduced in eq. (pq) . In addition, 
the identities (|104|) and (|101|) can be combined together to obtain: 



{da - 2dT,2f Ml = -2[m - m (m, r)] ds/n 



(106) 



By observing that the ratio {da — 2(i^2) /da is equal to (1 +r^), we can also express 
eq. ( |106| ) in terms of the pion decay constant F^, = —da and the quark condensate 
(ipip) = —ds/n. In this way we obtain: 



(1 + Fl Ml = 2{m-m) {^jji/j) 



(107) 



that is again eq. (^) derived before. Thus this equation, as its continuum analogous, 
is a direct consequence of the axial WI. 

We now show how the continuum current algebra relations and the partial conser- 
vation of the axial current are reproduced in the lattice theory at the leading order of 
the strong coupling and large expansion. Technically, the mechanism turns out to 



be exactly the same discussed in ref. |]14| for the weak coupling regime. 
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In the continuum theory, the axial WI, expressed in terms of renormahzed operators 
Oftix) and the renormahzed quark mass m^j has the form: 

d^MR,(^)niy)) = '^rnn {Phix)Phiy)) - ^ 5{x - y) (^(x)) (108) 

On the lattice, since the fermionic action is simply linear in the external source p{x), 
the derivatives of the partition function entering in the WI (p7| ) are proportional to the 
corresponding correlation functions. Thus, eq. (^) can be also written in the form: 

Vj^(Aj(a;)P^(?/)) = 2 (m + 4r) {P^{x)Pi{y)) - 8r (n^(x)P-''(y))- 

1 (109) 
n .y ^ \ ^' 

Formally, this differs from the continuum expression ( |108| ) for the presence of terms 
proportional to r, originated by the Wilson term in the lattice action. We now observe 
that these contributions can be expressed in terms of correlation functions of the two 
lattice operators V^A^{x) and P{x) respectively^. Explicitly, from eq. (pOf), we find: 



= (I-Za) V]f{Al{x)P^{y)) - 2m{P\x)P3{y)) 



;iio) 



where: 



Z^ = l-2^ = l + r2 (111) 

Therefore, by substituting eq. ( |110| ) in eq. ( |109| ), we obtain for the lattice WI the 
following expression: 

ZAVl%Al{x)P^{y)) = 2 (m - m) {P\x)P^{y)) - -5'^ 5,,y{S%x)) (112) 

Despite the strong coupling limit we are considering, this result has exactly the form 
( [108| ) required by continuum current algebra, provided we identify the following renor- 
mahzed operators and quark mass: 

A^^(x) = Z^Aj(x) , P'^{x)=ZpP\x) ^^^^^ 

S%{x) = ZpS^{x) , rriR = Zm {m - m) 



^In the weak coupling region this is also true, but for the presence of additional Schwinger terms. 
This follows from the fact that V^A^(a;) and P{x) are the only operators of dimension four or less 
with which the Wilson terms can mix 1141. 
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with Za given by eq. ( |111| ) and = 1/Zp. In particular, the vanishing of the renor- 
mahzed quark mass identifies the proper chiral hmit, in which the axial current 
becomes exactly conserved and the pseudoscalar meson mass vanishes. 

To verify the extent to which chiral symmetry is in fact recovered also in the strong 
coupling limit, we can study other examples of lattice WI's, like for instance the axial 
WI for the (ASP) correlation function. In this case, current algebra requires: 



V^"(4(x)5-'(i/)P'(^)) = 2(m-m) ^^{P''{x)S^ iy)P'iz))~ 



ZjnZf 

'z'a 



;il4) 



An explicit calculation of the correlation functions in the effective theory shows that 
again the WI has the form prescribed by eq. (|1 14|) . Za is still given by eq. (|111|) and 
Zm = ^/Zp. In addition, we also find the ratio: 

Zs , ^ 

/ = 1 (115) 

Zjp 

As a last example, we wish to discuss a vector WI. As far as the vector current is 
concerned, we know that, for any value of the gauge coupling, this current is conserved. 
Thus, its renormalization constant must be exactly equal to 1. This can be verified by 
considering, for instance, the vector WI for the (VPP) correlation function. In this 
case, in the limit of degenerate quark masses, the prescription of current algebra is: 

ZvV^^^V;:ix)P^{y)P\z)) =if'^^5,,y {P^{x)P\z)) + ^^^^^ 

+ tf^''^5.,z{P'{x)P^{y)) 

This is exactly the identity we obtain in the effective theory, where, as expected, we 
also find: 

Zv = 1 (117) 

To conclude this section, we wish to notice that above results provide a clear inter- 
pretation of eq. ( |107| ). Since the combination Z^Zp is equal to 1, the right-hand side of 



eq. ( |107| ) is renormalization invariant. In contrast, the pion decay constant P^- on the 
left-hand side of the equation must be renormalized, the proper renormalization con- 
stant being Za = 1 + r"^. Thus, once expressed in terms of all renormalized quantities, 
with Ft^fi = ZaFt, = Fq, eq. ( |107|) becomes: 



F^j,M^ = 2mR{^^lj)n (118) 
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that reproduces the well known continuum relation. 

Thus, all the results obtained in this section are consistent with the following picture: 
at the leading order in the strong coupling and large N expansion, the lattice theory 
reproduces correctly the relations of continuum current algebra. This allows one to 
define, even in the strong coupling limit, renormalized lattice operators that corresponds 
to their continuum counterpart. The lattice vector current, being conserved, does not 
require renormalization. In contrast, the axial current does. The resulting renormalized 
operators then respect all requirements of the PCAC hypothesis. The way in which this 



scenario sets up is exactly the one outlined in ref. |]T4| for the weak coupling regime of 
lattice QCD. 

6 The "renormalized" effective Lagrangian 

The discussion of the two previous sections has been mainly devoted to investigate the 
role of the Wilson terms (eq. (^2])) in the lattice effective Lagrangian. 

The original motivation, which led to the introduction of the Wilson term in the 
lattice action, was the removing of the doubler fermion species |T^. In the naively 
discretized version of the Dirac theory, these fermions would originate additional pseudo- 
Goldstone bosons. Clearly, were these particles present in the spectrum of the theory, 
they could not have been integrated out from the effective Lagrangian. Thus in this 
sense, by having considered only a multiplet of pseudo-Goldstone bosons, this primary 
effect of the Wilson term has been taken into account, form the very beginning, in our 
lattice effective theory. 

A well known consequence of removing the doubler species is the appearance on the 
lattice of the axial anomaly [|l^]. However, in the effective theory considered in this 



paper, we have not introduced external sources coupled to flavour- singlet anomalous 
axial currents. Thus, the Wilson terms in our effective action do not play any role in 
reproducing the correct anomaly. From this particular point of view, they could have 
been completely neglected. 

Beside that, most of the other effects of the Wilson term on the lattice are usually 
undesired, in that they contribute to further differentiate the lattice theory from the 
corresponding continuum one. By being formally vanishing in the continuum limit, 
the Wilson term in the original lattice action can only induce physical effects through 
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the loops of the quantum theory. In the weak couphng region, these effects can be 
cured by considering a proper set of renormahzation prescriptions, which allows one to 
relate the lattice quantities to the corresponding continuum counterparts. Examples 
of such prescriptions, already encountered in this paper, are the renormahzation of 
the quark mass and the multiplicative renormahzation of the bilinear quark operators. 
The several Wilson terms in the effective lattice Lagrangian enter to exactly reproduce 
these "undesired" effects of the Wilson term in the correlation functions of the lattice 
operators. 

In this context, it appears natural to consider instead the following approach: on 
one hand, one can simply discard the Wilson terms in the lattice effective theory, thus 
ending up with an effective Lagrangian with a reduced number of terms each of those 
has a direct correspondence in the continuum chiral Lagrangian. On the other hand, 
in comparing the predictions of this effective theory with the results of a numerical 
simulation, one must consider, in the latter, only the correlation functions of properly 
renormalized operators. 

In the following, we will show how this approach works in details by considering two 
specific examples. In doing that, we will take advantage of the fact that, as discussed in 
the previous sections, from the particular point of view of renormahzation, the lattice 
theory at the leading order in the strong coupling and large N limit shows basically 
the same features expected in the weak coupling region. Thus, we will be able to 
continue our discussion by dealing with the effective lattice Lagrangian derived in this 
approximation. 

The quark mass renormalization 

In the limit of vanishing bare quark mass and Wilson term, the generation of a quark 
mass from pure quantum corrections in the lattice theory is protected by chiral sym- 
metry. In contrast, in the presence of the Wilson term, a physical quark mass can arise 
even in the limit in which the mass term of the lattice action is set to zero. This is a 
well known consequence of the Wilson breaking of chiral symmetry. 

In this respect, the local external source x(x) of the lattice Lagrangian cannot be 
identified, directly, with the corresponding source of the continuum theory. In fact, 
when in the latter this source is removed, the pseudoscalar Goldstone bosons remain 
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exactly massless. On the lattice, instead, when we consider the limit x(^) = 0) cor- 
responding to m = — 4r, we still find a finite value of the pion mass (precisely, from 
eq. (ID, Ml = -16rV(l + r2)). 

In order to define an external source that is proportional, in the continuum limit, to 
the corresponding source of the continuum Lagrangian, we can perform in the effective 
action the following change of variable: 

-^(^x) = {m + Ar) ^{x) + 6x{x) (119) 

where m is the quantity defined in eq. (|105|) and the Wilson source ^(x) is introduced 
in eq. (|119|) to preserve chiral covariance. The new field Sx{x) is defined in such a way 
that, in the limit of vanishing external sources, it reduces to the renormalized quark 
mass: 

(5x)o = = m — m (120) 

For simplicity, we are only considering here the additive renormalization of the quark 
mass, thus neglecting the existence of an additional multiplicative renormalization. 

By construction, after this change of variable the whole pion mass in the effective 
theory is generated by the new source 5x- Precisely, as a consequence of the substitution 
( [119| ), a new term appears in the lattice effective action, of the form: 



-NC^{m + Ar) Y^Ti [i\x) U {x) + i{x)U\x) 
2 ^ 



:i2ii 



The contribution to the pion mass coming from this term is exactly canceled by an the 
opposite contribution coming from the Wilson terms. As a result, the pion mass in the 
effective theory is only generated by the term Sm of eq. (p4D, now containing the new 
source 5x- 

If in a specific case we were only interested in the calculation of the pion mass, it 
is clear that, after the change of variables (|119D, we could completely forget about all 



the Wilson terms in the effective action. In addition it is worth to notice that, but 
for the formal substitution x{x) ~^ ^xi^)) the mass term Sm in the effective action 
is left unchanged. Thus in this respect, we could have also considered, from the very 
beginning, a completely different point of view. That is, we could have discarded all 
the Wilson terms in the lattice effective action and simply interpreted the source x(^) 
as the "renormalized" source associated with the renormalized quark mass 5m. 
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Technically, there is a difference between the two particular points of view. This is 



due to the fact that, after the substitution ( |119| ), the Wilson sources a{x) and tt{x) are 
coupled in the action to different operators, namely and 11' (x) given by: 



S'(x) = S(a;) + (m + 4r) S{x) 
n'(x) = n(x) + (m + 4r) P{x) 



;i22) 



However, as far as we are not interested in calculating the correlation functions of these 
particular operators, the above difference becomes completely immaterial. 

Other interesting consequences of the change of variable defined by eq. ( p.l9| ), and 
in particular those related to the lattice chiral WI's, will be discussed at the end of this 
section. Before that, we want to consider another significative example. 



The axial current renormalization 

In the previous example we have shown that some of the Wilson terms in the effective 
Lagrangian can be removed by considering an external source which is directly coupled 
to a renormalized quantity, in that specific case the quark mass. Similarly, we can argue 
that other Wilson terms would be also removed by considering for instance a new axial 
source, a'^{x), directly coupled to the renormalized axial current, Aji^^^x) = ZaA^{x). 
Specifically, we would like to consider a change of variable leading to the introduction 
of the following set of external sources: 



W'^ix) 



- exp |- 
exp{- 



v'Jx) 



a^[x) 



(123) 



Since the vector current is conserved, one might look for the new vector and axial 
sources in the form: 



a^{x) = ZAa Jx) 



(124) 



with = given by eq. ( |1 1 1| ) in the strong coupling limit. 

This would imply a relation between the old and new and Z^ external sources, 
which takes the form of an infinite series: 



W^ix) 



r 
~2 



Wl{x) 



Z'^yx] 



W'^{x] 



+ . . 
+ . . . 



(125) 
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where the dots represent terms which are at least quadratic in the external sources v'^ 
and a'^. However, the changes of variables of eqs. ( |124| ) and (|125| ) cannot be implemented 
as such, not being consistent with chiral covariance. Requiring the validity of eq. ( |124| ) 
is too restrictive, since adding to the right-hand sides of that equation any term linear 
in the Wilson external sources, (y{x) and vr(a;), would not change the desired result. 
Indeed, even with this addition, the sources v'^{x) and a'^{x) would still be coupled, 
in the lattice action, to the currents V^{x) and ZaA^{x) respectively. Similarly, we 
could also add any other term which is at least quadratic in any of the defined external 
sources. 

Taking advantage of this freedom one can obtain chirally covariant relations between 



old and new sources, satisfying all the desired requirements, by modifying eq. ( 125 ) in 
the following manner: 



^2 



w,{x) = wlix) + - me{x)wi{x) + wi{x)i{x + m\x +n)- 



2e(x)Z'(x)et(x + /2) 



^ ■ ■ ■ (126) 



Z,{x) = Z'^ix) + - [^{x)^i{x)Z'^{x) + Z'^{x)i\x + ri)i{x + fi)- 



-2e{x)W'{x)ax + fl) 



+ 



In this expressions we have included on the right-hand side all the terms which are 
at most linear in the v'^, a'^, a and vr external sources, whereas the dots represent 
an infinite series of terms at least quadratic in these sources. The missing terms can 
be determined, order by order, by requiring the unitarity of the new fields VF^(x) and 
Z'^{x). By expanding eq. ( |126| ) in terms of vector, axial, scalar and pseudoscalar sources, 
we find, up to higher orders: 



^A^) = v'Jx) +0{f 



r 



a,ix) = (1 + r') a'^ix) - y V^7r(a;) + Oif) 



(127) 



which generalize eq. ( |124| ). 



It is worth to discuss an interesting consequence of the combined change of variables 
( [119| ) and ( |127|) . Once both these changes are performed, the terms in the effective 
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action which are hnear in the external sources can be still written in the form (^) and 
(pop, but now with new coefficients di given by: 



d, = -nNCm , dsi = NCm/F^ 

dp = ~NCm/FQ , da = —FQ , dv = 1 (128) 

da = 8nr^ [NC, - F^/ (1 + r^)] , d^i = 
d^i = , d^2 = 

In particular, because of the last two equalities {d.„i = d.„2 = 0), the operator: 

n'(a:) = n(x) + (m + 4r) P{x) + ^Vf^A^X^) (129) 

that is now coupled to the source '7r(x) in the lattice action, has always vanishing on-shell 
matrix elements, and can only contribute to the lattice correlation functions through 
localized contact terms. In some cases it does not contribute at all. For instance, the 
lattice chiral WI for the (AP) correlation function can be now directly derived in the 



form (|112|) , without any explicit additional contribution from the Wilson terms. 



The continuum-like effective Lagrangian 

The above discussion naturally leads to the following conclusion: since it is possible, 
both in the weak and strong coupling limit, to renormalize the lattice correlation func- 
tions in such a way that they satisfy all the relations of continuum current algebra, it 
must be also possible to describe the physics of this renormalized lattice theory in terms 
of a "renormalized" , continuum- like, effective Lagrangian. Specifically, the only sources 
of chiral symmetry breaking in this Lagrangian must be represented by the quark mass 
terms, whereas Wilson terms like those in eq. ( ^2]) should not enter at all. 

The two examples discussed in this section were aimed to support the above argu- 
ment. An evident feature of eq. (|128| ) is that, after the external sources renormalization, 
many of the coefficients related to the Wilson operators S(x) and 11 (x) do actually van- 
ish. In addition, the Wilson parameter r does not explicitly enter anymore in the 
expressions of the remaining coefficients, related to the lattice operators V^, ZaA^, S 
and P. In fact, one can easily show that after the change of variables (|119|) and (|126|) 
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many of the correlation functions we have considered so far are simply reproduced by 
the following effective action: 



S{U) = J2Tt \If^ {A,U)^ {A,U) - ^NCm {x^U + U^x 



(130) 



provided the external sources Wfj,{x), Z^{x) and x{^) entering in it are identified with 
the renormalized sources of eqs. ( |119| ) and ( |123| ). Equation ( |130| ) exactly corresponds, 
up to discretization effects, to the continuum chiral Lagrangian at order . 

On the other hand, it should be pointed out that a continuum-like effective action, 
like that of eq. (|130|) , cannot be obtained from the original one by simply renormalizing 



the external sources. For instance, the correlation functions of the operator A^^{x) 
obtained from eq. (|130|) are not equal to the corresponding correlation functions given 
by the lattice effective Lagrangian after the renormalization of the axial source. The 
reason is that the lattice the operator A^^{x) does not renormalize with renormalization 
constant and in fact it does not even renormalize in a simple multiplicative way. 
Only some of the effects of the Wilson terms in the lattice effective theory can be 
accounted for by a renormalization of the external sources, because these effects rather 
imply the renormalization of an infinite number of lattice operators. The important 
point here is the following: with a proper renormalization of the lattice operators, their 
correlation functions will be correctly reproduced, up to pure discretization effects, by 
a continuum-like effective lattice Lagrangian which does not contain the Wilson terms 
at all. The reason is that these correlation functions satisfy all the WI's predicted by 
continuum current algebra, and these identities can be only reproduced by a continuum- 
like effective Lagrangian. 

7 A numerical calculation of the lattice effective La- 
grangian in the scaling region 

In the region of intermediate couplings, which is relevant for numerical simulations of 
continuum QCD, the non-effective degrees of freedom of the lattice theory cannot be 
analytically integrated out. Therefore, in this case we would not be able to derive 
the exact form of the resulting effective theory. However, such an integration can be 
performed numerically, at least in an approximate way. 
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The idea consists in defining an effective lattice Lagrangian of the general form de- 
rived in this paper for the strong coupling and large limit. The basic fields entering 
in this Lagrangian are the effective field U{x) and the set of external sources Wfj_{x), 
Z^{x) and x{^)- Oii the basis of the results derived in this paper, we expect that the 
correlation functions of properly renormalized lattice operators can be reproduced by 
an effective lattice Lagrangian which does not necessitate Wilson terms. Furthermore, 
the external sources of this effective Lagrangian have a direct and well defined corre- 
spondence with the external sources entering in the continuum QCD chiral Lagrangian. 

The most general form of the effective lattice Lagrangian is dictated by local chi- 
ral invariance, lattice symmetries, parity and charge conjugation. The chiral paral- 
lel transport between two nearest neighbor lattice sites is performed by the left- and 
right-handed external sources Z^{x) and Wf^{x) respectively. They can enter in the 
Lagrangian in the form of chiral plaquettes or through the lattice covariant derivative 
A^U{x) defined in eq. (P7|). The local source x{^)y which reverses chirality, allows to 
introduce in the theory the effects of chiral symmetry breaking induced by the light 
quark mass term. 

The general strategy to calculate the effective lattice Lagrangian has been outlined 



in ref. [|I0[- The idea is to assume a sufficiently large set of effective couplings, with 
strength determined by unknown numerical coefficients. These coefficients can be then 
fixed through the matching of an overcomplete set of expectation values, computed 
both in the effective and the full theory. 

The calculation in the full theory, being highly non perturbative, requires the im- 
plementation of a numerical lattice simulation. An open question is whether or not 
the lattice effective Lagrangian, in the region of intermediate values of the lattice QCD 
coupling constant, can be treated in perturbation theory. This is possible in the strong 
coupling and large limit, where the coefficients of the higher order terms (0^, 0^, 
. . .) are suppressed by increasing powers of 1/A^. In addition, even for the actual value 
A^ = 3, these coefficients (see e.g. eq. (0)) are found to lie in the perturbative regime. 
Certainly, such a feature could be lost as one approaches the weak coupling region. 
Even if this should be the case, however, the difficulty could be more technical than 
conceptual, in the sense that also the effective theory could be treated by numerical 
techniques. 

A crucial observation is that, since the lattice effective Lagrangian contains explicitly 
the collective fields responsible for the long distance behaviour of the fundamental 
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lattice theory, then only short distance couplings are expected to play a relevant role in 
the effective theory. Therefore, in defining a criterion to truncate the infinite number 
of allowed couplings in the effective Lagrangian one can limit oneself to considering 
local or quasi-local interactions (e.g. only local, nearest and next-to-nearest neighbors 
couplings). For the same reason, the numerical simulation in the fundamental theory, 
which is necessary to perform the matching, should be feasible on lattices of moderate 
sizes, which would allow one to achieve a higher numerical accuracy. 

In the strong coupling limit, the coefficients of the effective Lagrangian are 
found to be explicit functions of the light quark mass m and the Wilson parameter r. 
Both these parameters enter in the coefficients through the vacuum expectation value 
Uq. Such a dependence persists even after the process of renormalization, but there 
are reasons to believe that it is an artifact of the strong coupling expansion. In the 
weak coupling limit, the matrix elements of the lattice renormalized operators must 
be independent on the value of r. Therefore, the coefficients of the renormalized effec- 
tive Lagrangian cannot depend on the Wilson parameter. Furthermore, there is clear 
phenomenological evidence that the effects of the light quark masses, in the physical 
amplitudes, are well reproduced in the chiral QCD Lagrangian by the terms contain- 
ing the external source x(a:), without any additional mass dependence entering in the 
coefficients. For these reasons, in defining an effective lattice Lagrangian for the weak 
coupling region one should assume the coefficients independent on m and r. 

Once one has been able to determine an expression for the lattice effective La- 
grangian which accurately reproduces the Green's functions of the fundamental lattice 
theory at moderate distances, one can compute its large volume limit and expand in 
increasing powers of the external momenta. The result of such an expansion, to any 
given order in p, corresponds to the continuum low-energy QCD effective theory. In 
such a way, the procedure discussed in this paper would allow a first principle theoretical 
calculation of the coefficients of the QCD chiral Lagrangian. 
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